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Martensitic  phase  transitions  with  surface  effects 


MARK  LUSK 

Division  of  Engineering  and  Applied  Science,  California  Institute  of  Technology,  Pasadena,  CA 
91125,  USA 

Abstract.  A  model  of  martensitic  phase  transitions  is  presented  that  is  capable  of  accounting 
for  a  variety  of  surface  effects  associated  with  the  localized  interaction  of  coexisting  phases  of 
a  material.  Such  phenomena  are  thought  to  play  a  critical  role  in  determining  the  size,  shape 
and  stability  of  nucleated  embryos  as  well  as  to  affect  the  conditions  under  which  nucleat'on 
occurs.  Attention  is  restricted  to  transitions  that  are  mechanically  induced,  and  the  model 
does  not  account  for  temperature  effects.  Materials  that  undergo  martensitic  phase  changes 
are  modeled  as  hyperelastic  in  both  the  bulk  and  the  interface.  The  characterization  of  such 
bodies  is  examined  in  detail,  and  a  representation  theorem  is  presented  for  isotropic,  hyperelastic 
interfaces  in  terms  of  two  scalar  invariants  associated  with  the  deformation  of  the  interface. 


1.  Introduction 

A  subject  of  continuing  interest  in  both  the  mechanics  and  materials  science 
communities  is  non-difFusive  phase  transitions  that  involve  two  solid  phases  of  a 
material  separated  by  a  sharp,  coherent  interface.  These  martensHic  processes 
are  most  common  in  Fe-C  systems,  but  axe  found  in  Ni-Ti,  Ag-Cd,  Au-Cd,  Cu- 
Zn,  Cu-Zn-Al,  Cu-Al-Ni,  In-Tl  and  other  alloys  as  well  [1-7].  In  addition,  some 
ceramic  materials  such  as  partially  stabilized  zirconia  exhibit  this  type  of  phase 
transition  [8-10] . 

In  recent  years  the  materials  science  community  has  developed  a  wide  vari¬ 
ety  of  products  whose  working  properties  rely  on  martensitic  phase  transitions, 
luese  include  shape-memory  materials  like  the  Nitinol  and  Tinel  alloys  [11].  An¬ 
other  technology  based  on  martensitic  phase  transitions  is  that  of  transformation- 
toughened  materials  [8-10]. 

Interest  in  martensitic  materials  has  motivated  a  variety  of  theoretical  in¬ 
vestigations  that  help  to  better  understand  these  phase  transitions  by  modeling 
them  from  a  continuum  physics  perspective  [12-34].  Here  phases  refer  to  disjoint 
domains  of  a  single  energy  functional  that  characterizes  the  material  of  interest. 


This  is  in  line  with  the  intuitive  picture  thivt  individual  phases  are  just  differing 
configurations  of  a  single  substance  and  provides  a  means  of  distinguishing  the 
deformation  of  a  given  phase  from  a  phase  transition.  The  classical  example  of 
this  type  of  material  description  is  that  for  the  Van  der  Waals  fluid  {35|.  From 
the  continuum  point  of  view,  an  interface  is  coherent  if  the  multiphase  configura¬ 
tion  can  be  described  as  a  continuous  deformation  of  some  homogeneous  reference 
state. 


Material  constitution  generally  imposes  severe  restrictions — phase  segrega¬ 
tion  requirements — on  the  ways  coherent  coexistence  can  be  accommodated.  Mul¬ 
tiphase  states  can  be  supported  only  under  certain  conditions,  and  the  sha])e  of 
each  phase  is  often  restricted  as  well  [34].  The  interface  between  phases  may 
accrete — that  is,  move  relative  to  the  underlying  material  so  as  to  effect  a  trans¬ 
fer  of  mass  from  one  phase  to  another.  However,  it  is  not  a  free  boundary  because 
of  the  phase  segregation  requirements,  and  certain  materials  may  call  for  supple¬ 
mentary  information  regarding  interface  position  [21-23,  26,  29,  30). 

Some  continuum  models  do  not  capture  certain  key  phenomena,  known  as 
surface  effects,  eissociated  with  phase  transitions.  These  eure  due  to  the  localized 
interaction  of  coexisting  material  phases.  This  interaction  is  particularly  relevant 
in  situations  where  the  interface  curvature  is  very  large — as  in  the  initial  formation 
of  a  new  phase  embryo  [41]. 

One  such  surface  effect  is  the  existence  of  a  jump  in  the  traction  exerted 
on  either  side  of  an  interface  within  a  statics  setting.  In  fluids,  for  example,  the 
pressure  inside  a  vapor  bubble  is  higher  than  that  of  the  surrounding  liquid,  but 
in  solids  it  is  reasonable  to  think  that  interfaces  may  resist  shear  as  well,  since 
this  is  a  characteristic  of  the  bulk  phases  on  either  side  of  the  interface. 

A  second  important  effect  attributable  to  interface  properties  is  the  occur¬ 
rence  of  supercritical  phenomena  [36-40].  Conditions  capable  of  initiating  a  phase 
change  with  planar  interface  geometry  cannot  produce  a  phase  embryo  with  a 
highly  curved  boundary.  For  example,  a  phase  change  might  not  occur  until 
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the  body  is  superheated  above  its  normal  transition  temperature,  Likeunse,  for 
mechanically  induced  reactions,  a  body  may  require  straining  beyond  the  level 
known  to  induce  phase  changes  when  inhomogeneities  edlow  the  fomiation  of  low 
curvature  interfaces. 

The  symmetry  properties  associated  with  an  interface  influence  the  shape  of 
a  phase  embryo,  and  this  is  especially  true  when  the  nucleus  is  so  small  that  bulk 
effects  do  not  dominate  [41,  42],  The  surface  properties  of  an  interface  can  thus 
play  an  important  role  in  the  determination  of  embryo  shape. 

A  final  surface  effect  pertains  to  the  stability  of  newly  formed  nuclei  which 
may  depend  upon  the  embryo  size  as  well  as  on  conditions  at  the  system  bound¬ 
ary.  This  behavior,  and  often  the  existence  of  such  equilibria  themselves,  are 
attributable  to  localized  properties  at  phase  interfaces  [38-40]. 

Continuum  treatments  of  phase  transitions  may  be  modified  to  account  for 
these  important  surface  effects.  Interface  properties  are  modeled  by  fields  defined 
on  the  interface  itself,  which  are  collectively  referred  to  eis  interfacial  structure. 
Here  the  interface  is  usually  treated  as  a  surface.  Balance  principles  are  postulated 
that  restrict  the  way  in  which  bulk  and  surface  fields  may  interact,  and  materials 
are  characterized  by  constitutive  response  functions.  An  historical  survey  of  the 
evolution  of  such  models  is  provided  by  LuSK  [45]. 

For  simplicity,  attention  is  restricted  to  martensitic  transitions  that  are  me¬ 
chanically  induced,  and  the  development  presented  here  does  not  account  for  tem¬ 
perature  effects.  The  theory  is  a  modified  version  of  that  developed  by  Gurtin 
and  Struthers  [46],  and  offers  an  alternate  foundation  for  studying  surface  ef¬ 
fects  in  martensitic  phase  transitions.  Materials  that  undergo  such  changes  are 
modeled  as  hyperelastic  in  both  the  bulk  and  the  interface.  The  characterization 
of  such  bodies  is  examined  in  detail,  and  a  representation  theorem  is  presented 
for  describing  the  interfaces  of  isotropic,  hyperelastic  bodies  in  terms  of  two  fun¬ 
damental  scalar  invariants  associated  with  the  deformation  of  the  interface.  An 
application  of  the  theory  is  discussed  in  Lusk  [47]. 
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2.  Preliminaries 


Within  a  purely  mechanical  theory  of  ph^lse  transitions,  phase  interfaces  are  asso¬ 
ciated  with  strain  discontinuities.  For  clarity  of  presentation,  the  theory  presented 
here  assumes  there  exists  at  most  one  such  interface  in  a  given  body. 


2.1  Notation  and  definitions 


E^' 

Unit(V^) 

Lin(V^,V'''0 

Lin'^^(V'''',V'‘^) 

Lin+(V^’,W^) 

Sym(V^’) 

Sym+(V^) 

Unim(V^,i/^) 

Orth(V^,W^) 

Orth+(V^,fY^) 


set  of  all  A'-tuples  of  real  numbers 

iV-dirnensional  inner  product  spaces 

•V-dimensional  Euclidean  space 

set  of  all  unit  length  elements  of 

set  of  all  linear  maps  of  V‘^’  — » 

set  of  all  non-singular  elements  of  Lin(  V'^'  ,  V''^) 

set  of  all  elements  of  Lin(V^,i/^)  with  positive  determiiifint 

set  of  all  symmetric  elements  of  Lin{  ) 

set  of  all  positive  definite  elements  of  Sym(  ) 

set  of  all  elements  of  Lin^{V^ M^)  with  unit  determinant 

set  of  all  orthogonal  elements  of  Lin(  ) 

set  of  all  orthogonal  elements  of  Lin+(V^,W'’' ) 


The  inner  product  of  a,  b  G  is  denoted  by  <  a,  b  >  or  a  •  b.  The  outer 
product  of  c  €  and  d  €  is  denoted  by  c  0  d  €  Lin(  V'^ ,  ).  The 

vector  product  of  f,  g  €  is  denoted  by  (f  A  g)  6  V^.  The  inner  product  on 
Lin(V'^,  V^)  is  defined  by 

A  •  B  :=  Trace(AB^)  =  (A,B)  V  A,  B  G  Lin(V^,  (2.1) 

so  that  all  of  the  sets  of  linear  transformations  defined  above  are  endowed  with 
an  inner  product  space  structure. 

In  general,  scalar  fields  are  denoted  by  lower-case  letters,  vector  fields  by 
lower-case  bold  type,  and  tensor  fields  by  bold-type,  upper-case  letters.  A  super- 
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script  denotes  partial  differentiation  with  respect  to  the  second  argument  of  a 
function. 


Surface  S 

Surface  normal  n(x) 

Tangent  space  n-^(x)  of  S  at  x 
Surface  boundary  dS 

Tangential  vector  field 
Tangential  tensor  field 
Boundary  bi-normal,  m(x) 

Projection  map  P  of 

Inclusion  map  I  of  5 


a  subset  of  locally  diffeomorphic  to  E‘. 

either  of  the  two  unit  vectors  perpendicular 

to  a  surface  at  the  point  x. 

best  linear  approximation  to  S  at  x. 

the  closed  curve  that  delineates  the  edge  of 

a  surface. 

field  on  S  with  values  in  n‘^(x). 
field  on  S  with  values  in  Lin(n-^{x).  n-^(x)) . 
the  outward  pointing  unit  vector  perpendic¬ 
ular  to  dS  at  x  but  within  the  tangent 
space  of  S  at  that  point. 

tensor  field  on  S  with  values  in  Lin{JE'^.  n-^(x)) 
defined  by  P(x)a  =  a  -  (a  ■  n(x))n(x) 

Vx  €  5,  Va€ 

tensor  field  on  5  with  values  in  Lin(n^(x),  iE^) 
defined  by  I(x)a  =  a  Vx  G  5,  Va  G  n'^(x). 
P(x)I(x)  =  I5,  the  tangent  space  identity 
map  at  x,  while  I(x)P(x)  =  1  -  n(x)  ©  n(x). 
Also,  I  =  P^. 


2.2  Bulk  and  surface  gradients 

Let  g  :  E^ .  Then  for  every  x  in  the  gradient  of  g  at  x,  denoted  by 

Vg(x),  is  the  unique  element  of  Lin(iE^,iE^)  such  that 

lim  (^?(x  -h  hk)  -  {Vg{x))hk^  =  g{x)  Vk  G  E^ .  (2.2) 

Now  consider  a  scalar  field  g  defined  on  S  instead  of  JE^.  Let  (f>  parameterize 
S — that  is  4>  :  E^  —*S  diffeomorphically.  Then  there  exists  /  :  E^  —*  E  given  by 


(2.3) 


the  composition  /  :=  go<l>  which  allows  the  gradient  of  g  to  be  defined  as 

Vg(x):= 

It  is  straightforward  to  show  that  this  definition  of  the  gradient  of  a  scalar  field 
is  independent  of  the  parameterization  0  of  S.  The  idea  is  extendable  to  vector 
and  tensor  fields  on  .S  as  well.  For  emphasis,  denote  by  the  surface  gradient 
of  a  field  on  S-  Note  that  if  w,  b.  and  A  are  scalar,  vector,  and  tensor  fields  on 
S,  respectively,  then 

V^u>(x)  €  n-'(x) 

V^b(x)  €  Lin(n^(x),iE;'‘)  (2.4) 

V.:  A(x)  €  Lin(n‘^(x).Lin(2E'^,  lE^)). 

2.S  Smooth  fields  and  surface  divergence  ' 

Let  M  :  S  M  he  a.  map  such  that  ^  u;  is  defined  on  S.  Then  w  is  smooth  on  S 
and  is  this  is  denoted  by  u?  €  C^(S)  or  sometimes  by  w  :  S  ^  IR.  The  extension 
to  vector  and  tensor  fields  is  completely  analogous.  If  b  is  a  smooth,  vector  field 
on  S  then  ^-b  :=  Trace(P^  b)  on  .S'  is  the  surface  divergence  of  b.  Analogously, 
for  C  a  smooth,  tensor  field  on  S  the  surface  divergence  of  C  is  given  by 

(V^-C)k  =  VkeiE;^  (2.5) 

The  tangential  tensor  field  L  :=  — n  provides  curvature  information  in¬ 
trinsic  to  the  surface  and  is  therefore  called  the  curvature  tensor.  It  can  be  shown 
that  L  is  symmetric  and  that  the  scalar  field  H  ;=  ■|Trace(L)  gives  the  surface 
mean  curvature.  A  corollary  to  the  Stokes  theorem  that  is  useful  for  surface 
considerations  is  the 


Surface  Divergence  Theorem. 

Let  f  be  a  smooth,  vector  field  on  S.  Then 


J  V,  ■  fdA 

s 


j  {  mdL- 

dS 


n)dA. 


(2.6) 
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2-4  Surface  motion 


Consider  a  one-parameter  fsunily  of  surfaces  «S(f)  defined  over  some  open  time 
interval  T.  Let 

:={(x,t)|x  e  5(<),  teT] 

(2.7) 

(55)t  :={(x,t)|x  €  teT} 

define  the  surface  and  surface-boundary  trajectories.  Let  <j>  :  U  C  JE^  xT  —*Sr  be 
a  one-parameter  family  of  surface  par ameterizat ions.  Denote  by  the  restriction 
of  <f>  to  dU  X  T,  where  dU  is  the  boundary  of  U .  If  there  exists  a  4>  such  that  (p 
is  defined  on  U  x  T  and  (f>d  is  defined  on  dU  x  T,  then  the  surface  is  said  to  be 
smoothly  propagating.  It  can  be  shown  that  while  the  interface  and  edge  velocities 


V(x,t)  :=0(^-‘(x,t),t) 
Va(x,t)  :=^e{4>'\x,t),t) 


(2.8) 


depend  on  the  choice  of  parameterization  <i>  the  normal  and  intrinsic  edge  speeds, 


V„  :=V  •  n  on  Sr 


(2.9) 


Vm  :=V  .  m  on  dSy, 

are  independent  of  parameterization.  V„n  is  the  normal  velocity  of  S.  Vniin  is 
the  intrinsic,  tangential  edge  velocity  of  S.  On  the  surface  boimdary,  the  sum  of 
these  two  velocity  fields  is  referred  to  as  the  intrinsic  edge  velocity  of  S. 

Consider  a  normal  trajectory  of  Sr  through  x  6  5{<o),  €  T.  This  is  a 

smooth,  t-parameterized  curve  in  lE^  described  by  a  set  of  position  vectors  z  such 


that 


z{t)£Sit),  P(z(<))i(f)  =  0,  VteT.  (2.10) 


The  normal  trajectory  is  used  to  define  a  time  derivative  following  S  by 

?(x,0  :=  (2.11) 

where  z  gives  the  normal  trajectory  of  Sr  through  x  6  *S(f).  It  may  be  shown 
that 


n=-V,F„. 


(2.12) 
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This  derivative  also  finds  application  in  the 
Surface  Transport  Theorem. 

Let  fii  be  a  smooth  scalar  field  on  a  smoothly  propagating  surface  S.  Then 

^  J  gdA=  J  l^-2HgVr,]dA+  j  gV^dL  Vt  €  T.  (2.13) 

s(t)  S(t)  dS(t) 

2.5  Latticed  bodies  and  interfaces 

A  collection  B  of  elements  p  is  a  body  if  there  is  a  set  Q{B)  o{  continuous  bijections 
X  such  that: 

(i)  the  gradient  of  the  induced  map  between  any  two  elements  has  positive 
determinant  where  it  is  defined; 

(ii)  Q.{B)  contains  all  continuous  bijections  of  its  elements. 

A  lattice  C{B)  of  a  body  B  is  a  set  of  elements  of  such  that: 

(i)  the  gradient  of  the  induced  map  between  any  two  elements  is  smooth; 

(ii)  C{B)  contains  all  smooth  bijections  (diffeomorphisms)  of  its  elemems. 

A  latticed  body  is  any  body  together  with  one  of  its  lattices.  The  view  is 
taken  that  a  body  must  be  provided  with  a  lattice — a  preferred  group  of  reference 
configurations — in  order  to  identify  phase  interfaces  kinematically. 

A  two-phase  deformation  with  respect  to  ^(B)  is  a  configuration  A*  G  ^(B) 
such  that  the  induced  map  between  fC  and  X  is  arbitr8irily  smooth  everywhere 
except  possibly  on  a  phase  interface  S.  A  two-phase  motion  X^.  is  a  one-parameter 
family  of  two-phase  deformations  A:'(#)  6  Q{B)  with  t  ranging  over  an  open  interval 
T.  It  is  admissible  if  there  exists  a  reference  configuration  K.  G  n(B)  such  that 

(i)  each  element  of  Xj-  is  a  two-phase  deformation  with  respect  to  K.\ 

(ii)  Sr  is  a  smoothly  propagating  surface. 

For  a  latticed  body  to  undergo  an  admissible,  two-phase  deformation  it  must  be 
admissible  with  respect  to  one,  and  hence  cill,  of  the  elements  of  its  lattice. 
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Consider  the  two-phase  deformation  shown  in  Fig.  1  with  induced  map  y 
and  bulk  deformation  gradient  F  :=  Vy  defined  where  y  is  smooth.  Then  the 
surfaice  in  the  deformed  configuration  of  the  body, 

•5  :=  {y|y  =  y(x),  x€5},  (2.14) 

is  referred  to  as  the  deformed  interface  with  {ii,  m}  and  P  its  surface  orientations 
and  projection  map.  The  restriction  of  y  to  denoted  by  y,  is  the  invertible 
deformation  map  of  the  surface  with  an  cissociated  surface  deformation  gradient 
F  •=  ^y  defined  on  S.  Note  that  F(x)  6  Lin+ (n-^(x),  n-^(y )),  y  =  y(x), 
even  though  F  is  not  defined  on  S.  It  is  often  convenient  to  use  the  linear 
tremsformation  representing  the  surface  gradient  of  y: 

JF:=  V^y  =  F+I  =  F-I,  (2.15) 


where 


Thus, 


F-  =  jim[F(xl/n(x))]  Vx  e  S. 
F  =  PF,  F-  =FP-h(F-)n(gn, 


and  it  is  easy  to  show  that 


F  =  iF. 


(2.16) 


(2.17) 


(2.18) 


Note  that  F  e  Lin'"^(n^,F='). 

The  orientations  of  the  undeformed  and  deformed  interfaces  may  be  related 
using  the  determinants  of  the  deformation  gradients 


J  :=Det(F) 
j  :=Det(F) 


Det(F) 

iF-^n|| 


(2.19) 


where  F  may  be  either  F"'"  or  F  in  the  second  equality.  The  relations  are 


J 


n, 


m  = 


m, 


(2.20) 
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where  F  may  be  either  F"^  or  F“  and  where  h  is  the  determinant  of  the  gradient 
of  the  restriction  of  y  to  dS. 


Define  the  material  and  spatial  velocity  fields,  v  and  v  respectively,  as 

v{x,t)  :=^(x,t) 


v(y,t)  :=::v(y-*(y,f),t). 

If  V  is  a  surface  velocity  of  Sf,  then  the  induced  interface  velocity 

V  :=  (Fl)V  +  v^ 


(2.21) 


(2.22) 


is  a  surface  velocity  for  Sr — the  interface  trajectory  associated  with  deforming 
interface  S.  As  (2.22)  indicates,  the  definition  of  V  does  not  depend  on  the  side 


of  S  used  and  in  fact 


V  =  FcV  +  Vc  Vc€[0,l] 


(2.23) 


where 


Fc  :=cF-^  +  (l-c)F- 
Vc  +  (1  —  c)v" . 


(2.24) 


Let  and  K.2  be  lattice  elements  related  by  the  linear  transformation  G  : 
Ki{S)—*K2{B).  If  Vi  is  an  interface  velocity  for  the  motion  with  respect  to  ACj, 
then  V2  =  GVi  is  an  interface  velocity  for  the  motion  with  respect  to  K.2  with 
V]  and  V2  inducing  the  same  interface  trajectory  Sr- 


2.6  Galilean  objectivity 

Consider  the  admissible,  two-phase  motions  Xr  and  of  a  latticed  body  {B,C} 
with 

X*it)  =  QX(t)-hdt  +  e  VteT, 

and 

Q  €  Orth+(iE;^  d,  e  6  E\ 


-11- 

Denote  by  Hit)  :=  X{t;B)  the  region  occupied  by  B  for  all  t  in  T,  with  Sr  the 
deformed  interface  trajectory  partitioning  'R{t)  over  T.  Define  the  bulk  trajectory 
as 

:=  {(y>0|y  €  teT}. 

Suppose  that 

4>  •  (^T  \  ^ 

A:  {nr\Sr)-*Un{E^,E^) 

^  :  Sr  —*  E 

A  :  (y,0  6  Sr  •->  A(y,t)  €  Lin(h-^(y,  t),  JE^) 

are  fields  associated  with  {B,£}  during  motion  AV  and  that 

A*;  (7l;  \5;)-^Lin(^^JE:^) 

:  s;-^R 

A*  :  (y,0  G  5;  A-{y,<)  €  Lin(iV(y, t),  iE^") 

are  the  analogous  fields  of  {5,£}  during  motion  The  fields  associated  with 
{B,C}  axe  said  to  exhibit  Galilean  objectivity  if 

=  (t>{yj) 

V(y,<)  EEr\Sr 

A*(z*(y,0,<)  =  QA(y,0Q" 

and 

0*(z*(y,t),f)  =  ^(y,t) 

V(y,t)G5T, 

A*(z*(y,<),t)  =QA(y,t)Q^ 

where 

z-ly,*)  ;=  Qy  +  dt  +  e  V(y,t)  € 
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3.  Balance  postulates 

3.1  Primitive  fields 

Consider  a  latticed  body  (IS,P}  and  an  admissible  two-phase  motion  .\V.  and 
let  T^r  and  Sr  be  the  bulk  and  interface  trajectories  discussed  in  section  2.0. 
Stipulate  that  for  every  such  motion  the  following  smooth  fields  exist: 


Bulk  mass  density 

i>  ■ 

{TZr  \Sr)-JRr 

Bulk  energy  density 

ir  : 

{nr\$r)--IR 

Cauchy  stress 

T  : 

{itr\Sr)^lui{E\E^) 

Interfacial  energy  density 

IV  : 

Sr^m 

Deformational  surface  stress 

t  , 

(yj)  €  Sr  >—  t{y.t)  6  Uu{n^{y,t).  E^) 

Accretive  surface  stress 

C  : 

(y.O  6  5t  t(y.f)  €  Lin(n-^{y. 0. 

3.S  Balance  postulates 

Let  Vr  be  the  bulk  trajectory  associated  with  an  arbitrary  subbody  of  B  under 
the  previous  motion  AV  as  shown  in  Fig.  2.  Denote  by  Qj  the  partitioning 
subinterface  trajectory  associated  with  Pr.  The  boundary  of  Q  is  denoted  by 
dQ.  Let  0  be  the  outward  unit  normal  to  dP — the  boundary  of  P.  Interface 
normal  h  and  bi-normal  til  are  as  defined  in  section  2.  The  following  postulates 
are  imposed  on  all  such  subregions  and  motions.  The  theory  presented  here  is  a 
modified  version  of  that  developed  by  Gurtin  and  Stru'I  HERS  [’6]. 


Mass  Conservation 
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Momenta  Balances 


£ 

dt 


j  pv  dV  —  j  Ti/dA-i-  j  Tiii 


dL 


V(t) 


SQit) 


dt 


J  y  A  pv  dV  =  j  y  ATi^dA+  J  y  A  Tin  dL 


T(() 


dQ(t) 


Mechanical  Dissipation  Imbalance 


£  +  A'  <  P, 


where 


E{t-V{t))  ;=  J  WdV+  J  wdA 

r{t)  C(() 

K{LV(t))  :=  j  ipv-vrfV 


Vit) 

•{t;V{t)):=  J  Tu-vdA+  J  Tm-ydL+  j  Cm-(V-v+)dI 

dV{t)  aQ{t)  dQ(t) 

+  J  Cm  - (V  -  v-)d£. 

dQ(t) 


Parameterization  Invariance 

The  power  P  expended  at  the  boundary  of  P  is  independent  of  the 
parameterization  of  Q  for  all  times  in  T. 

Objectivity 

W ^  w,  T,  T,  and  C  exhibit  Galilean  objectivity  as  defined  in  section  2. 

By  Cauchy’s  theorem,  the  momenta  balances  imply  that  T  is  symmetric.  A 
curvilinear  triangle  version  of  the  theorem  applied  on  the  interface  shows  that  the 
deformational  surface  stress  T  is  symmetric  as  well,  and  is  therefore  representable 
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as  a  tangential  tensor  field  as  defined  in  section  2.1.  Henceforth  T  is  taJken  to  be 
tangential.  Note  that  symmetry  of  T  and  T  together  with  conservation  of  linear 
momentum  implies  the  angular  momentum  balance. 

The  presence  of  two  distinct  interfacial  stresses  taxes  physical  intuition,  and 
this  issue  is  addressed  at  several  points  in  the  ensuing  development.  The  defor- 
mational  stress  T  expends  power  in  association  with  the  absolute  movement  of 
the  interface.  Even  during  a  two-phase  motion  for  which  no  new  material  changes 
phase,  work  is  performed  by  this  stress.  In  contrast,  the  accretive  stress  C  ex¬ 
pends  power  in  association  with  the  relative  motion  of  the  interface  with  respect 
to  the  underlying  material — a  process  referred  to  as  accretion.  In  the  limiting 
case,  not  covered  by  this  theory,  of  a  phase  transition  occurring  with  no  material 
deformation  the  accretive  stress  may  still  expend  power.  This  stress  is  the  Eule- 
rian  counterpeirt  to  that  presented  by  Gurtin  and  Struthers  [46]  and  is  related 
to  the  capillary  tractions  of  Leo  and  Sekerka  [48],  Cahn  and  Hoffman  [49]  and 
Hoffman  and  Cahn  [50].  Its  physical  significance  for  hyperelastic  materials  is 
elucidated  in  sections  4  and  5. 

3. 3  Referential  formulation 

The  derivation  of  local  field  equations  is  facilitated  by  re-expressing  ail  primitive 
fields  as  quantities  defined  on  an  arbitraxy  element  of  the  body  lattice,  as  shown 
in  Fig.  2.  To  this  end,  define  the  following  fields  on  such  a  reference  configuration: 

Referential  mass  density  p  :—Jp 

Referential  bulk  energy  W:=JW 

Piola-Kirchhoff  stress  S  :=JTF"^ 

Referential  interfacial  energy  w  :=jw 

Referential  deformational  surface  Stress  S  :=;TtF-^,  jf  =  PSF^ 

Referential  accretive  surface  stress  C  •.=2j((F^))CF"^, 
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where  F,  F,  j,  and  J  are  as  defined  in  section  2.5  and 


((<^{x)))  :=  lim  i{^(x  + /in(x))  +  <ji(x  - /in(x))  }  Vx  €  (3.1) 

h'\0  Z 


Using  these  fields,  the  linear  momentum  balance  and  the  mechanical  dissi¬ 
pation  imbalance  can  be  expressed  as 


^  j  pvdV  J  SudA+  J  SmdL 

T  dV  dQ(i) 


(3.2) 


and 

£ 

dt 


J  W  dV  +  ^  J  ^  dt  J  ^ 


Q(<) 


(3.3) 


<  J  Su-vdA  +  J  Sm  •  [FcV  +  Vc]  dL  +  J  Cm  -Vdi, 

dV  dQ{t)  dQ{t) 

with  Fc  and  defined  by  (2.24). 

The  respective  localizations  of  the  linear  momentum  balance  away  from  and 
local  to  the  interface  are 

V  •  S  =pv 

Ipv|U„  +  [SnJ  =  -  ^  •  S. 

Here 


(3.4) 


[<;^(x)J  :=  lim{^(x  +  hn(x))  -  ^(x  - /in(x))}  Vx  G  5.  (3.5) 

h'^O 

Analogous  localizations  for  the  dissipation  imbalance  are  made  after  considering 
the  field  restrictions  imposed  by  the  principle  of  parameterization  invariance. 

3-4  Implications  of  parameterization  invariance 

Use  of  the  surface  divergence  theorem  gives 

y  Cm  •  V  dl  =  y  V,  •  (C^m  ■  V;„m)dA, 


dQ{t) 


Q(0 


(3.6) 
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and  application  of  the  product  rule  then  yields 


j  Cm-y  dL=  J  [-C-LV;-C^n-ii-f(V^-C)nV;]d^+  J  CmmV^dL. 


Now  consider  the  second  term  in  the  referential  dissipation  inequality, 


J  Sm  •  [FcV  +  Vc]  dL.  (3.S) 

ac(0 

By  the  surface  divergence  theorem  this  is  equal  to 

J  y,-lS^F,nV„]dA+  J  V,  [S^Vc]dA+  J  Sm  •  [F,mV^]dL 

C(«)  Q{i)  dQ(t) 

=  J  Sm-(F,m)V^dI-  J  [{FS)  ■  L  +  {%  ■  S)  ■  iF,n)]V„  dA  ^3 

aQ(t)  Q(t) 

+  J  |(SP).f',-(S^F,n).il+(V,.S)-v,|d^. 

Q(f) 

The  mechanical  dissipation  imbalance  is  now  expressible  in  a  form  amenable  to 
application  of  the 


Invariance  Lemma.  Let 


e  ;  (^T  \  ^t)  — ^ 

/  :  Sr 

g  :  St -^F, 


and  suppose  that 


JedV+  J  fdA+  J  9V,ndL<0  'iV  CTI,  'it  £7. 


Then  g  =  0  on  Sr  ■ 


This  is  proved  in  [46].  The  result  restricts  the  primitive  fields  by  the  parameteri¬ 
zation  invariance  postulate,  which  comes  in  the  form  of  the 
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Tension-Energy  Theorem  [46]. 


F^S-f  PC  =  u>l^. 


(3.10) 


Equations  (3.3)  and  (3.7)-(3.10}  then  combine  to  give  a  form  of  the  mechan¬ 
ical  dissipation  imbalance  that  is  consistent  with  the  parameterization  invariance 
postulate: 


jwdV+j^  j  wdA+jJ^p’^ydV 


dt 


V  BQd)  V 

<  j  js  l'c-ac  n-bV^  S  v,|d.4-l-  J 


Su  •  V  dA 


Q{t) 


dv 


+  J  wV^dL+  J  [(V^-S)-(Fcn)-H(V^-C)  n-2^u;]V;dA, 


dQit) 


QU) 


where 


ac  :=  (S^Fc  +  C^)n. 


(3.11) 

(3.12) 


3.5  Referential  localizations  and  interface  driving  traction 

Application  of  the  surface  transport  theorem  and  linear  momentum  balance  to 
(3.11)  yields  the 

Reduced  Dissipation  Imbalance 


£ 

dt 


JwdV<  Js  FdV+  J  [(V,-S)-((Fn))-h(V,-C)-n]V„dA 


e(«) 


-  y  [iS  -  S  •  f’c  -h  ac  •  in] . 


(3.13) 


Q(t) 


Localized  away  from  the  interface,  this  implies  that 


S  .  F  -  IF  >  0. 


(3.14) 
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The  local  implication  at  the  interface  is 

fV„  -  [u-  -  S  -  +  ae  •  n]  >  0,  (3.15) 


where 

/  ~  I'l-l  -  {(S))  ■  [F]  +  (V,  .  S)  ■  ((Fn))  +  (V,  .  C)  n  (3.16) 


is  the  interface  driving  traction.  Equation  (3.15)  was  first  derived  by  Knowles 
[15]  within  a  setting  devoid  of  surface  fields.  Gurtin  and  Struthers  [46]  derived 
(3.15)  in  its  present  form,  but  their  work  differs  from  the  present  development  in 
that,  in  [46],  a  more  extensive  set  of  postulates  leads  to  a  local  balance  principle 
with  the  form  of  (3.16). 

Equations  (3.4),  (3.10),  (3.14)  and  (3.15)  along  with  the  angular  momentum 
localizations, 

SF^  =  FS’’ 


SF^  =  FS^, 


(3.17) 


comprise  the  local  field  equations  for  the  mechanical  theory  of  dissipative  accre 
tion. 
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4.  Hyperelastic  materials  with  interfaces 
4.1  Material  constitution 


The  mechanical  theory  summarized  by  (3.4),  (3.10)  and  (3.14)-(3.17)  is  incom¬ 
plete,  since  a  means  of  characterizing  the  bulk  and  interface  of  a  given  materiEil 
has  yet  to  be  provided.  This  is  done  by  establishing  constitutive  relationships 
between  the  primitive  fields  and  the  process  kinematics.  Though  the  form  of 
such  relations  may  be  quite  general,  attention  for  the  remeunder  of  this  work  is 
restricted  to  hyperelastic  materials  with  interfaces.  For  such  materials,  there  exist 
energy  response  functions  W  and  u'  such  that 


W  =  W^(F) 
S  =  aK''(F) 


(4.1) 


w  =  wiFc,  n) 

SP  =  6,ii>(Fc,  n) 


(4.2) 


ac  —  5jU’(Fc>  n), 


with 

W  :Un+{E\E^)^lR 


(4.3) 

w{-,n):Un+{E\E^)^lR  Vn  G  Unit(JE:^), 
and  where  di(f>  refers  to  the  partial  derivative  of  <f>  with  respect  to  its  i**  argument. 
The  symbol  d<i>  refers  to  the  derivative  of  a  fimction  <t>  of  a  single  tirgument.  Here 
and  in  the  following  development,  P  and  I  are  treated  as  elements  of  \j\n[E^ ,  it‘‘  ) 
and  Lin(n-^,iF^),  respectively,  defined  by 


Pb  =  b~(b-n)n,  I  =  P^,  VbGiS^  (4.4) 


The  projection  and  inclusion  maps  P,  I  associated  with  the  deformed  surface 
normal  n  are  treated  analogously. 

For  such  materials,  (3.14)  is  strictly  satisfied,  implying  that  the  bulk  material 
does  not  dissipate  energy.  The  interface  is  still  dissipative,  however,  with  (3.15) 
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reducing  to 

fV„  >  0,  (4.5) 

where  the  driving  traction  /  is  given  by  (3.16).  Accretion  is  therefore  the  only 
way  that  hyperelastic  materials  can  dissipate  energy,  as  is  clear  from  both  (4.5) 
and  the  fact  that 

f:  +  A'-P  =  -  /  fVndA  (4.6) 

G(0 

for  such  materials.  Here  E,  K,  and  P  are  as  defined  in  the  mechanical  dissipation 
imbedance. 

The  interface  response  functions  of  (4.1)  and  (4.2)  can  be  expressed  in  a  more 
useful  form  using  the  following  lemma  from  [46]: 


Lemma  1.  Consider  a  function  w  defined  by 
tj;(iF,  f,  n)  =w{FP  +  f  0  n,  n) 

Vn  €  Unit(^='),  VfeiE:’,  VF  6  Lin'"®(n^,  iEJ^). 

Then 

diw{lF,  f,  n)  —diw{]FP  +  f  (g)  n,  n)I 
diw{F,  f,  n)  =5,it;(FP  +  f  <8)  n,  n)n. 


(4.7) 


(4.8) 


Lemma  1  and  (4.2)2  imply  that 

S  =  5,ii;(F,  f,  n)  (4.9) 


and  that 

a,iD(F,f,n)  =  0.  (4.10) 

The  surface  fields  are  therefore  functions  only  of  n  and  F,  and  the  surface  energy 
of  (4.2)  could  just  as  well  have  depended  upon  F"*"  or  F";  preference  has  not  been 
given  to  the  bulk  on  either  side  of  the  interface  in  describing  the  interfacial  fields. 

A  final  step  in  the  material  characterization  is  to  express  the  referential  accre¬ 
tive  stress  C  in  terms  of  the  surface  energy  response  function  w.  This  is  facilitated 
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by  introducing  the  partial  derivative  of  w  with  respect  to  n  following  the 

interface,  identified  as  the  unique  element  of  n'^  such  that 

D„w{F,  n)  •  a  =  ±u,(Fq{k(0),  11),^,^)) 

d0  \  ^  J  (4.11) 

Vk  ;  JR^UnitCiE*),  k(0)  =  n,  k{0)  =  a, 

where 

Q  :  Unit(^:^)  X  Unit(iE:^)  ^  Orth+(iE^ iE’) 
is  the  bilinear  map  defined  by 

Q(e,  f )e  =  f,  Q(e.f)g  =  g  V{e,f)  €  Unit(iE^)  x  Unit(iE^), 
with  g  •  e  =  g  •  f  =  0.  Then,  as  shown  in  [46], 

Dr,w{F,i,n)  =  d,w{FF  +  f®n,n)  +  [drw{FJ,n)Yi-F’^d^w[F,i,ny  (4.12) 
Since 

C  =  IPC  +  n  (8)  C^'n,  (4.13) 

(3.10)  and  (3.12)  and  the  fact  that  PI  =  combine  to  give 

C  =  wl  -  IF^'S  +  n  (ac  -  S^Fcii) .  (4.14) 

Equations  (4.2),  (4.7),  (4.9)  and  (4.12)  then  yield 

C  =  w{F,  f,  n)I  -  \F'^d,w{F,  f,  n)  -  n  g)  £>„ie(F,  f,  n).  (4.15) 


Equations  (4.9),  (4.10)  anvl  (4.15)  admit  the  following  characterization  of 
hyperelastic  materials  with  interfaces; 


w  =  iy(F) 

S  =  dW{F), 


(4.16) 


and 


w  =  w{F,  n) 

S  =  dtW{F,  n) 

C  =  w{F,n)l  -  lF^diw(F,n)  -n®  DnU>(F, n). 


(4.17) 
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Note  that  the  response  of  the  system  is  independent  of  the  scalar  c  introduced  in 
(3.3).  The  tangential  component  of  C  refers  to  PC,  and  only  the  first  two  terms 
in  the  above  expression  for  C  are  represented  in  its  tangential  component  since 
Pn  =  0.  The  last  term  in  (4.17)3  is  referred  to  as  the  normal  component  of  the 
accretive  stress  and  is  equivalent  to  n  (g)  C^n. 

4-2  A  kinetic  relation  for  interfaces 


The  mechanical  theory  of  hyperelastic  materials  summarized  by  (3.4),  (3.17), 
(4.5)  and  (4.16)-(4.17)  is  still  incomplete.  As  shown  by  Ericksen  [12]  and 
Abeyaratne  and  Knowles  [21-23,  26]  for  the  case  of  no  surface  fields,  there 
are  settings  for  which  the  location  of  the  interface  Ccinnot  be  determined  uniquely 
for  prescribed  boundary  data.  Abeyaratne  and  Knowles  suggest  that  this  may  be 
interpreted  as  a  constitutive  deficiency.  Motivated  by  the  internal  variable  work 
of  Rice  [43,  44]  and  the  form  of  (4.5),  they  offer  a  simple  remedy  in  the  form  of 
a  constitutive  relation  between  the  driving  traction  /  and  the  accretive  velocity 
of  the  interface  V„: 


v;  =  v{f), 

v{f)f>o  v/e2R. 


(4.18) 


Equation  (4.18)i  is  referred  to  as  a  kinetic  relation  with  the  inequality  restricting 
V  imposed  by  the  dissipation  imbeilance  (4.5).  Supplementing  the  system  postu¬ 
lates  and  material  constitution  with  such  a  relation  has  been  shown,  in  certain 
settings,  to  resolve  the  uniqueness  problem  [23,  26,  28-30].  Though  developed 
without  considering  surface  fields,  the  rationale  for  adopting  a  kinetic  relation  is 
still  valid  within  settings  that  include  interfacial  structure.  Such  a  constitutive 
remedy  is  therefore  adopted  in  the  present  work.  It  should  be  noted,  however,  that 
Gurtin  and  Struthers  [46]  have  derived  an  equation  with  the  form  of  (4.18)i 
from  a  more  elaborate  set  of  balance  postulates. 
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^.5  Mechanical  equivalence 

Consider  two  lattice  elements  ATj  and  /C2  with  respect  to  which  a  body  is  hyper- 
elastic,  that  are  related  by  linear  transformation  H,  where  K2  =  HoATi.  These 
two  reference  configurations  re  said  to  be  mechanically  equivalent  if  first, 

p,(x,t)  =  p2(Hx,t)  Vx  €  7^l(^)\5l(t),  VteT  (4.19) 

and,  second,  they  are  characterized  by  the  same  response  functions: 

IT  2  ~  I'i  ],  UI2  =  •  (4.20) 


By  virtue  of  the  construction  of  referential  fields  given  in  section  3.3,  the 
following  formulae  relate  mass  density  and  energy  fields  independent  of  any  con¬ 
stitutive  considerations: 


P2(x,0  =  Det(H)pi(Hx,t) 

1 


IT^2(x,t)  = 


Det(H) 


M^i(Hx,f) 


VxG7^2\5l(t),  VteT 


(4.21) 


and 


U)2(X,t)  = 


iH^n(x,i 


-i(;i(Hx,t)  Vx  6  5i(t),  VteT. 


(4.22) 


Det(H) 

Equations  {4.19)-(4.22)  collectively  imply  that  two  reference  configurations  are 
mechanically  equivalent  if  and  only  if 


Det(H )  =  1 

Ty^(F)  =  W"i(HF)  VFe  Lin+(JE^.2E:=’),  Vn  e  Unit(iE:^). 

(4.23) 

Define  the  material  symmetry  group  of  a  latticed  body  with  respect  to 
lattice  element  fC  to  be  the  set  of  all  lattice  elements  that  are  mechanically 
equivalent  to  fC.  If  Qk:  —  Unim(iE'^,  jE*^),  the  material  is  called  a  fluid.  If 
=  Orth+(iE^,  iE^),  the  material  is  called  a  solid. 
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The  kinetic  relation  does  not  play  a  role  in  determining  the  mechanical  equiv¬ 
alence  of  two  lattice  elements.  An  appropriate  chsinge  of  configuration  formula 
can,  however,  be  established.  Let  J„  eind  ju  be  the  bulk  and  surface  determinants 
associated  wfith  the  linear  transformation  H.  Then 

(4.24) 

3h 

and 

/i  =  J,j2  (4.25) 

so  that 

V2(4-)=—Viif)  V/6  2R.  (4.2G) 

Jh 

4-4  The  role  of  accretive  stress 

Recall  that  the  local  power  expenditure  associated  with  C  is  given  by  Cm  •  V. 
For  hyperelastic  materials,  then,  (4.17)3  implies  that 


Cm-V=  {ra-  [(F^a,ii;)m]  [(J9»-m]r„.  (4.27) 

The  first  term  represents  power  expended  in  extending  the  interface  boundary  via 
tangential  accretion.  The  second  term  represents  a  power  expenditure  accompa¬ 
nying  changes  in  interfacial  orientation  in  association  with  normal  accretion.  The 
accretive  traction  has  a  physically  meaningful  relation  to  the  symmetry  exhibited 
by  the  interface  and  is  embodied  in  the  following 

Theorem  On  Fluid-Surfaces. 

(A)  For  fluids: 

(i)  w  =  a,  3.  constant; 

(ii)  C  =  0;  and 

(iii)  T  =  ^1^. 
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(B)  For  a  hyperelastic  bodies  that  do  not  support  an  accretive  stress: 

(i)  ti!  =  CT,  a  constant;  and 

(ii)  T  =  als. 

(C)  For  a  hypereleistic  material  with  constant  interfacial  energy  per  unit  area: 

(i)  T  =  d-1^;  and 

(ii)  C  =  0. 

(A)  is  due  to  GuRTlN  and  Struthers  [46],  while  (B)  and  (C)  are  new. 

Proo/.(B)  For  C  =  0,  the  tension-energy  theorem  of  section  3.4  implies  that 

F^S  =  u;l;,.  (4.28) 

In  terms  of  T,  this  is  equivalent  to 

=iF^ifF-"‘  =;F^TF-’'. 


Therefore, 

T  =  — F"’^!  F^  =  iDl 

j  ^ 

Equation  (ii)  is  therefore  proved  once  (i)  is  established.  It  may  be  shown  that 
(4.28)  is  equivalent  to 

S  ==  tolF-^.  (4.29) 

But  for  hyperelastic  materials,  (4.17)2  gives 

S  =  diw{F,n).  (4.30) 

Now  let 

cr(F, n)  :=  —w{]Fyn) 

so  that 

diw(lF,n)  =  ;TF"^cr(F, n)  +  jd\a{F,  n).  (4. 31) 


Here  use  has  been  made  of  the  identity, 

aJ(F,n)  =  jiF-’'  (4.32) 

from  [46),  with  j{lF,  ti)  =  j.  Equations  (4.29)-(4.31 )  then  imply  that 


u-{F,n)IF-^  =  jcr(F,n)iF-^  +;c),d(F.n). 

Therefore, 

n)  =  0. 

C  =  0  also  implies  that  C^n  =  0  which,  by  (4.17)3,  gives  that 

0  =  Dnw(lF,n)  -  [Z)nj(F,n)]cr(F.n)  +  J(F,n)[Dnd(F.n)]. 


where 


j(F.u) 


Det(F) 

IP-Tn!i 


from  (2.19)2.  But  following  the  definition  of  D„{  )  given  by  (4.11), 


T>n;(F,n) 


L 

d0 


d 

dl3 


d 

d0 


Det[FQ(k(^),n)]  ] 

||[FQ(k(^),n)]-^k(/?)|ij  ^=0 
Det(F)  \ 

l|F-^Q(k(/?),n)k(/?)j|  j  ^=0 
Det(F)  ] 

IF-^nll 


=  0. 


Since 

D„]iF,  n)  =  DnjiFP  -I-  f  0  n,n) 
by  construction,  (4.34)  and  (4.35)  imply  that 


DnuiF,  n)  =  0. 


(4.33) 


(4.34) 


(4.35) 


Use  of  (4.7),  (4.8),  and  (4.12)  with  a  in  place  of  w  then  gives  that  0  must  be  a 
constant-valued  function.  Thus,  w  —  a,  a  constant,  □ 


Proof. (C)  w  =  a  is  equivalent  to  w  —  jcr  so  that  S  =  jdlF~^ .  Expression  of 
this  result  in  terms  of  T  yields  T  =  and  also  implies  that 

IF^S  =  jal  =  tel. 

From  (4.30)  and  the  fact  that  w  =  ja, 

Dnw{F,  n)  =  0. 

Substitution  of  this  result  and  (4.33)  into  (4.17)3  yields 

C  =  U7l  —  tel~n®0  =  0.  □ 

Part  (A)  indicates  that  fluids  cannot  support  an  accretive  stress,  though 
the  converse  of  this  does  not  hold  in  general.  Parts  (B)  and  (C),  on  the  other 
hand,  do  not  depend  upon  the  mechanical  symmetry  of  the  body.  They  imply 
that,  independent  of  the  constitution  of  the  bulk,  an  interface  exhibits  a  fluid 
nature — that  is,  surface  energy  and  surface  tension  characterized  by  the  same 
scalar  constant — if  and  only  if  C  =  0.  Thus,  the  accretive  stress  accounts  for 
special  properties  of  interfaces  exibiting  a  solid  nature  that  are  not  associated 
with  their  fluid-like  counterparts. 
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5.  Isotropic,  hyperelastic  materials  with  interfaces 

A  standard  result  in  continuum  mechanics  asserts  that  hyperelastic  materials  are 
isotropic  with  respect  to  a  given  reference  configuration  if  and  only  if  the  bulk 
energy  can  be  represented  by  a  function  of  the  three  fundamental  scalar  invariants 
of  the  right  Cauchy-Green  tensor  C  :=  F'^F  or  its  square  root  U.  These  scalar 
invariants  are  defined  as 

/i(U)  ;=  Trace(U) 

/2(U)  :=  \  [(Trace(U))^  -  Trace(U2)]  VU  €  Sym+(iE:^),  (5.1) 

/3(U)  :=  Det(U) 

so  that,  in  the  absence  of  surface  fields,  a  hyperelastic  material  is  isotropic  if  and 
only  if  it  is  characterized  by 

W^Wiluh^h) 

s  =  {diW  -b  /ia2l^)FU-^  -  (a2lF)F  -b  {hd3W)F-^  (5.2) 

T  =  ^{di\\r  +  iid2W)y  -  ~(a2^)v''  -b  (as  1^)1, 

h  •‘3 

where 

Ik  :=  /fc(U)  =  /fc(V),  U  :=  and  V  ;=  VfF^.  (5.3) 

Material  isotropy  also  imposes  restrictions  on  the  interfacial  response  func¬ 
tions,  and  the  elucidation  of  these  restrictions  is  central  to  this  work.  One  such 
restriction  is  on  the  accretive  surface  stress  C  and  is  a  direct  result  of  the  following 

Proposition  1.  Material  isotropy  implies  that  DnW  =  0. 

Proof.  By  (4.23),  isotropy  implies  that 

u;(F,  n)  =  u;(FQ,Q’'n) 

(5.4) 

VQ  €  Orth+(JE^.E:3)^  VF  €  Lin+(^E;^  iF^),  Vn  €  Unit(.E^). 


Recalling  (4.11),  let 


q  =  Q(k(^),n) 
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for  some 

k:  k(0)  =  n. 


Then, 

Q^{k(0),ti)n  =  k{0)  (5.5) 


so  that  (5.4)  implies  that 

tIi(^FQ(k(/3),n),k(/3)^  =t5(F,n)  V^.  (5.6) 

Therefore, 

^tf.('Fq(k(/3),.i),k{^A  =  ^a(F,n)  =  0. 

Thus,  by  the  definition  of  DnU>  given  by  (4.11),  DnW  =  DnW  =  0.  0 

Prop.  1  implies  that  the  accretive  stress  of  isotropic,  hyperelastic  materials 
does  not  have  a  normal  component.  This  is  clear  from  (4.17).  In  light  of  this 
result  euid  using  (4.27),  it  may  be  concluded  that  any  work  performed  by  C  to 
reorient  the  interface  must  be  due  to  material  anisotropy  since  the  last  term  in 
(4.27)  vanishes  for  isotropic  bodies.  Leo  and  Sekerka  [48]  conclude  that  such  a 
work  term  should  be  present  for  anisotropic,  hyperelastic  materials,  and  the  idea 
is  also  in  line  with  the  earlier  work  of  Cahn  and  Hoffman  [49]  and  Hoffman 
and  Cahn  [50],  who  adopted  a  capillary  traction  vector  to  account  for  anisotropic 
effects  in  non- deformable  media.  In  fact,  consideration  of  such  a  work  term  can 
be  traced  back  to  Herring  [41,  42]. 

Attention  is  next  turned  to  the  representation  of  interfacied  fields  in  terms 
of  the  scalar  invariants  of  the  surface  deformation  gradient.  An  isotropic,  scalar 
function  4>  is  characterized  by  the  property  that 


<^(Q^UQ,n)  =  <i(U,n) 

VU  G  Sym‘*'(n-^),  VQ  6  Orth+(n'^,  n"^),  Vn  6  Unit(iE^^). 


(5.7) 


Such  functions  admit  an  equivalent  representation  in  terms  of  the  two  fun¬ 
damental  scalar  invariants  of  their  first  argument — a  property  expressed  in  the 
following 
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Proposition  2.  If  is  8U1  isotropic,  scalar  function,  then  there  exists  a  function, 
<i>  :i  y.  Unit(iE'^)  — +  JR  such  that 


where 


and 


<i.(U,n)  =0(t{U),j(U),n) 

(5.8) 

VU€Sym+(n^),  Vn  €  Unit(,E;^), 

1:=  {(Ci,C2)|(Ci,C2)€iRi,  (4C2-C?)<o| 

(5.9) 

i{U)  :=  Trace(U),  j{ij)  ;=  Det(U). 

(5.10) 

This  proposition  is  proved  with  the  help  of  the  following  two  lemmas. 

Lemma  2.  (Ci)C2)  ^  coincide,  respectively,  with  i(A)  and  j(A)  for  some 
A  6  Sym‘^(n-^)  if  and  only  if  ((1,(2)  €  J,  with  1  as  defined  in  (5.9). 

Proof.  Suppose  A  €  Sym'^(n'^).  Then  A  has  two,  real,  positive  eigenvalues 
X  satisfying  the  characteristic  equation 


p(A)  =  +  i( A)A  -  j(A)  =  0. 

Application  of  the  quadratic  formula  then  reveals  that  (i(A),j{A))  €  J.  Con¬ 
versely,  suppose  (C11C2)  G  T.  Then 

—X^  +  Ci-^  ~  C2  =  0 


has  two  positive  roots  Ai,A2  and  Ci  =  Ai  -f  X2  and  ^2  =  AiA2.  Construct  A  € 
Sym'''(  )  such  that  it  heis  a  component  representation  in  its  principal,  rectilinear, 

Cartesian  coordinate  basis  of 

Ai  0  ■ 

0  A2 

□ 


[A] 


Clearly,  2(A)  =  and  j(A)  =  (2- 
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Lemma  3.  Let  A,B  €  Sym'*'(n'^ ).  Then  t{A)  =  i(B),  j{A)  =  j(B)  if  and  only 
if  there  exists  a  Q  €  Orth+(n-^,  n-^).  such  that  B  =  Q’^AQ. 

Proof.  Suppose  B  =  Q^AQ,  Q  6  Orth+{n-^,  n-^).  Then 

2(B)  =  (B,1.)  =  (Q^AQ,1.> 

=  (A,QQ")  =  (A,1,)  =  z(A). 

Also, 

j{B)  =  j(Q"  AQ)  =  j{Q^)j{A)j{Q)  =  j(A). 

Conversely,  suppose  2(A)  =  2(B)  and  j{A)  =  J'(B).  Then  the  eigenvalues  of  A 
and  B  coincide  since  they  are  both  given  by  the  characteristic  polynomial 

-X^  +  t\-j  =  0. 

Let  {61,62}  and  {65,62}  be  the  associated  principal  baises  for  A  and  B,  respec¬ 
tively.  Then  there  exists  a  Q  6  Orth+(n'^, n-*-),  such  that  Qe^  =  60,  a  =  1,2. 
Thus  , 

Be',  =  =  A,Q^e,  =  Q^Ae,  =  Q’' AQe'„  a  =  1, 2. 

implying  that  B  =  Q^AQ.  □ 

Suppose  there  exists  a  function  ;  J  x  Unit(.E®)  — »  2R  such  that 

</>(U,n)  =  ^(2(U),j(U),n)  VU€Sym+(n^),  Vn  e  Unit(iB^). 

Then  by  Lemma  3, 

.^(U.n)  =  .^(.(Q^UC!),j(Q’‘UQ),n) 

=  «(Q''t;q,n)  VQ€Orlh+(n"-,  !!•"). 

Therefore,  0  is  an  isotropic,  scalar  function.  Conversely,  if  (p  is  an  isotropic, 
scalar  function  then  by  Lemmas  2  and  3  there  exists  a  single-valued  function, 
4>  :  i  X  UnitfiE^)  iR,  such  that 

<^(U,n)  =  ^(2(U),j(U),n)  VU  €  Sym+Cn-^),  Vn  e  Unit(iE:''). 
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Prop.  2  is  therefore  established.  It  is  will  now  be  used  to  show  that  the 
interfaces  of  isotropic,  hyperelastic  materials  may  be  cliaracterized  by  two  scalar 
invariants  of  the  interface  deformation  as  given  by  the 


l.N'TERF.ACE  Representation  Theorem. 


Hyperelastic  materials  whose  bulk  is  characterized  by  (5.2)  are  isotropic  if  and 
only  if  the  interfacial  energt  can  be  represented  by  a  function  w  :  1  —*  IR,  such 


that 


u'(F,n)  =te(t(U),;(U)) 

VF€  Vn  e  Unit(F^), 


(5.11 ) 


where 

U  :=  =  Vf^. 


(5.12) 


Proof.  That  the  representation  given  by  (5.11)  is  possible  for  all  isotropic 
materials  is  proved  with  the  help  of  two  lemmas.  The  first  establishes  that  the 
interfacial  energy  of  isotropic,  hyperelastic  materials  may  be  characterized  by  an 
isotropic,  scalar  function  while  the  second  lemma  edlows  the  orientation  depen¬ 
dence  of  the  surface  energj'  response  function  to  be  dropped. 


Lemma  4.  For  isotropic,  hyperel2istic  materials  the  referentizd,  interfacial  energy 
can  be  represented  by  an  isotropic,  scalar  function 

ii>(-,n):  Sym'*'(n'^ )  — ►  iR,  Vn  6  Unit(iF^). 


Proof  Galilean  objectivity  implies  that 


u-(QFI,n)  =  u;(FI,n)  VQ  €  Orth+(iE’^ F^), 

Vn  €  Unit(F*), 


VF  €  Un^{E\E^), 

(5.13) 


with  I  the  inclusion  map  eissociated  with  ii'^. 

By  the  polar  decomposition  theorem,  there  exist  unique 


U  6  Sym'*‘(n-^),  Q  €  0rth4.(n'^, ), 
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such  that  F  =  QU  for  every  F  6  Lin+(n'^,n^). 

Consider  a  particular  linear  transformation 

Q  =  IQ^P  +  n0n. 

Note  that  this  Q  is  in  Orth+(iF^, since 

QQ^  =  (IQ^P  +  n  0  h)(iQP  +  h  ©  n) 

=  IQ^PIQP  4-  IQ^Ph  0  n  +  (n  ©  h)IQP  +  (n  @  n)(h  ©  n) 
=  IQ^  1-  QP  +  n  ©  n 
=  IQ^QP  +  n©n  =  IP  +  n®n  =  l. 

For  this  choice  of  Q, 

QFI  =  QiF  =  QiQU  =  (IQ"^P  +  n  ©  h)iQU 

=  IQ^^PIQU  +  (n  ©  n)lQU  =  IQ^^l^QU  =  lU 
Thus,  Galilean  objectivity  implies  that 

ui(FI,n)  =  iD(IU,n)  VF  €  Lm+(iE:^ iF^),  Vn  €  Unit{iE:^), 

with 

U  =  VF^JF  = 

and  where  F  and  F  are  as  defined  as  in  section  2.5. 

Now  suppose  that  Q  :=  IQP  +  n  ©  n  with  Q  €  Orth+(n-^, n-^). 
Q^Q  =  1  and  Q  =  PQI.  From  (5.13)  and  (5.14), 

u;(FI,n)  =  tI;(IU,n)  =  u;(QIU,n). 

But  for  the  Q  chosen, 

QIU  =  (IQP  +  n  ©  n)IU  =  IQPIU  =  IQU. 

(5.15)  and  (5.16)  together  imply  that 

u;(IQU,n)  =  ie(IU,n)  Vn  €  Unit(JE;*),  VU  €  Sym+(n^) 


(5.14) 

Then 

(5.15) 

(5.16) 


VQ  €  Orth+(n-^,n‘^). 


(5.17) 
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Equation  (5.17)  holds  for  all  hyperelzistic  bodies,  but  material  isotropy  im¬ 
poses  the  additional  restricition  that 


t7(FQ,Q^n)  =  u>(F,n)  VF  €  Lin+(^^ Vn  G  Unit(iF^), 

VQ  G  Orth+{iE:\iE^). 

Choose  Q  :=  IQP  n  0  n,  as  in  the  previous  claim,  and  note  that 


(5.1S) 


Q^n  =  (IQ^P-|-n0n)n=  n. 

Thus,  isotropy  implies  that  t7(FQ,n)  =  tZ>(F,  n)  for  such  a  Q  which,  by  the 
construction  of  in  (4.7),  gives  that 


txi(FI,n)  =  iy(FQI,  n). 


(5.19) 


FQI  =  F(IQP  +  n  0  n)I  =  FIQPI  =  FIQ. 


Equation  (5.19)  can  thus  be  written  as 


wiF,n)  =  wiFQ,n), 


(5.20) 


and  this  holds  for  every  F  G  Lin^^(n-^, iE^),  n  G  Unit^ C^),  Q  €  Orth+(n-^,  n-^). 
Since  lU  G  Lin^'^(n'‘',iE^),  (5.20)  implies  that 


ii;(IU,n)  =  tZ>(IUQ,n)  Vn  G  Unit(iE^),  VU  G  Sym+(n-^), 

VQ  €  Orth+(n'^,  n"^). 


(5.21) 


Equations  (5.17)  and  (5.21)  and  the  fact  that  the  inclusion  map  I  is  deter¬ 
mined  solely  by  the  unit  normal  n  admits  construction  of  an  isotropic,  scalar 


function 


;(-,n)  :  Sym+(n^) iR,  Vn  G  Unit(iE®), 


such  that  17(1X7,  n)  =  t7(U,  n). 


Lemma  5.  d3w{i,j,n)  =  0  V {i,  j,  n}  G  J  x  Unit(iE'^). 
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Proof.  Let  F*  =  FQ,  n*  =  Q^n,  and  I*  =  Q’TQ,  with  F  €  Lm+{E\E^), 
n  €  \Jmt{E^),  Q  €  OTth.^{E^,E^),  and  I  the  inclusion  map  for  n"^.  Let 

U*  =  v/(F*I*)^(F*I*) ,  U  =  v/(Fir(Fr). 

Then  U  and  U*  have  the  same  determinant  since 

;'(U*)  =  j  (P‘(F-^FT )  =  i(Q’-PF^FIQ)  =  /(U), 

where  both  determinants  must  be  positive.  U  and  U*  must  also  have  the  same 
trace.  This  is  true  since 

t(U*2)  =  (P'F'^F*!*,!^)  =  (Q’'PF^FIQ,1-) 

=  (PF^FI,l5)  =  t(U2), 

and 

t(U)=  (^(U2)  +  2;=(U)]‘'^ 

Material  isotropy  implies  that  u>(F,n)  =  ie(F*,n*),  and  in  terms  of  w  this 
means  that 

w(i(U),j(U),n)  =  ^(.(U-).j(U*).n-). 

But  U  and  U*  have  the  same  trace  and  determinant  so  this  equation  is  equivalent 
to 

n)  =  iD(t,j,n*)  Vn,n*  6  Unit(iC’).  □ 

Because  of  this  result,  the  orientation  dependence  of  ti;  may  be  dropped  and  it; 
considered  a  function  ii>  ;  I -^R. 

By  (4.23),  a  sufficient  condition  to  guarantee  isotropy  for  materi2ils  w'hose 
bulk  is  characterized  by  (5.2)  is  that 

w(FQ,Q’'n)  =  iI.(F,n)  VF  €  Lin+(,E^  JEJ’),  Vn  €  Unit(iE:=*), 

VQ  e  Orth+liFME:^). 

In  terms  of  w  this  is  equivalent  to 

iZt(i(U),j(U))  =  u>(i(U’),;(U*))  VF€Lin+{iE;\2E;"),  Vn  €  Unit(iE;"), 

VQ€Orth+(iE:^£:'), 
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where  U  and  U*  are  as  defined  in  Lemma  5.  But  this  equality  is  guaranteed  since 
U  and  U*  have  the  szime  trace  and  determinmt  by  Lemma  5.  □ 

The  interface  representation  theorem  is  therefore  established.  In  conjunction 
with  its  coimterpart  from  the  elasticity  theory  of  bulk  materials,  this  theorem 
implies  that  a  necessary  and  sufficient  condition  for  a  hyperelastic  material  to  be 
isotropic  is  that  (5.2)i  and  (5.11)  be  satisfied.  Note  that  the  interfaces  of  such 
bodies  are  completely  characterized  b.  two  scalar  invariants  associated  with  the 
deformation  of  the  interfaces  without  regard  for  interfacial  orientation. 

Equation  (4.17)2  gives  the  referential  deformational  surface  stress  S  as  a 
function  of  iu  so  the  above  result  can  be  used  to  express  S  in  terms  of  w.  In  the 
following  discussion,  a  pr  .  ndicates  the  directional  derivative  of  a  function. 
Also,  it  is  convenient  t'  r  .i.oduce  the  two-dimensional  right  Cauchy-Green  tensor 
:=  U^.  The  directional  derivative  of  w  with  respect  to  its  first  argument  is 
then  given  by 

(S,iB)  =  tu'(F,n)B  =[ait^(t,j)t'(U)-ha2u;(i,j);'(U)]$'(Q)r'(F)iB 

VJB  €  Lin"'®(n-",F^), 

where 

$(Q)  ;=  VQ  €  Sym+(n^) 

r(F):=F^F  VF€  Lin^^(n'^,F^). 

A  standard  mathematical  result  (Ciarlet  [51]  10-11)  is  that 

i'(U)D=t(D)  VU,D  e  Sym+(n^) 

/(U)D  =i(U)i(U-*D)  VU,D  €  Sym+(n^). 

Also, 

r(F)F  =  +  VF,B  €  Lin'"*(n-^,  F^), 

$'(Q)D  =  Q,D€Sym+(n^). 

(5  25}  implies  that 

$'(Q)[r'(F}iB]  =  iu-'(B^F-f  F^F), 


(5.22) 


(5.23) 


(5.24) 


(5.25) 


(5.26) 
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and  this  result  can  be  used  with  (5.24)i  to  obtain 

+  F^B)] 

-  {F\J~\B). 

Likewise,  (5.24)2  and  (5.26)  imply  that 

;'(U)$'(Q)r'(iF)iB  =  ij(U)(U-=,B’'iF+  F’'B> 

=  j(U)(FU-‘,iB>. 

(5.27)  and  (5.2S)  can  then  be  applied  to  (5.22)  to  obtain 

(S,®>  =d,uii,j){Fij-\B)  +  d2w{iJ)j{ij){FV-\B) 
VB  e  Lin'"^(n^,iE:^). 

Therefore, 

S  =  diw{iJ)F\J-'  ■hd2w{iJ)jiiJ)F\J-^ 

and 

t  =  -aiui(i,;)PFU-‘F^  4-a2U.(i,;)PFU-*F^. 

J 


(5.27) 


(5.28) 


(5.29) 


(5.30) 


(5.31) 


These  representations  can  be  simplified  by  using  the  left  polar  decomposition  of 
the  surface  deformation  gradient.  Given  F  =  QU  =  VQ,  then  V  =  FU“‘F’'. 
Therefore, 

PFU-‘F’'=V.  (5.32) 

A  simplification  is  also  possible  by  noting  that 


PFU-*F"^  =  1_. 


(5.33) 


Thus,  the  deformational  stresses  of  isotropic,  hyperelastic  interfaces  can  be  rep¬ 
resented  as 

S  =  i[diw{i,j)F\J~'-  +  j{\J)d2w{i,j)F~'^  (5.34) 

and 

T  =  -diw(i,j)V  +  d2w(i,j)l^. 

J 

A  representation  for  the  accretive  stress  is  then  obtained  using  (4.17)3. 


(5.35) 
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6.  Summary 


A  model  for  martensitic  phase  transitions  has  been  presented  that  is  intended  to 
capture  localized  interfacial  effects.  A  set  of  physically  reasonable  postulates  lead 
to  field  equations  and  jump  conditions.  The  notion  of  a  hyperelastic  material 
was  extended  to  include  bodies  that  support  structured  interfaces,  and  an  inter¬ 
face  representation  theorem  was  derived  for  isotropic,  hyperelastic  bodies.  The 
following  is  a  summary  of  the  interfacial  characterization  of  these  materials: 


w  =  w{i,j) 

s  =  i{(a,uOFu-'  +  oa»F-"'} 

T  =  -f  (52w)l5 

C  =  I{-(5,ti))U  -f  {w  ~  jd2w)ls], 

where 

t:=z(LJ)  =  z(V),  j:=;{U)  =  ;(V) 

U  :=  \/1^,  V  :=  V¥¥^. 


(C.l) 


Thus,  isotropic,  hypcrelastic  interfaces  are  completely  characterized  once  the  de¬ 
pendence  of  the  surface  energy  on  the  invariants  of  the  interfacial  deformation 
is  established.  Moreover,  a  necessary  and  sufficient  condition  for  a  hyperelastic 
material  to  be  isotropic  is  that  (5.2)i  and  (6.1)i  be  satisfied. 

An  application  of  this  theory  is  examined  in  Lusk  [47],  where  three  different 
interfacial  constitutions  are  examined  in  a  probiemi  concerning  nucleation  and 
growth. 
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